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Abstract. We considcr unitary cvolution of finitc bipartite quantum systems and 
study time dependence of purity for initial cat states - coherent superpositions of 
Gaussian wavc-packcts. Wc derivo cxplicit formula for purity in systems with nonzero 
time averaged coupling, a typical situation for systems where an uncoupled part of thc 
Hamiltonian is Liouvillc integrable. Previous analytical studies have been limited to 
harmònic oscillator systems but with our theory we are able to derive analytical results 
for general integrable systems. Two regimes are clearly identified, at short times purity 
decays due to decoherence whereas at longer times it decays because of relaxation. 
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1. Introduction 

Quantum mechanics is a linear theory and as such a superposition of solutions is also an 
admissible solution. Appearance of such superpositions at the quantum level and at the 
same time their absence from the macroscopic classical world has troubled scientists 
from the very beginning, an example is the famous "Schrodinger Cat" p. Their 
absence is usually explained in terms of decoherence due to the coupling with external 
degrees of freedom, see j2] for a review. In recent years the process of decoherence 
has actually been experiment ally measured 0]. In the present paper we will derive 
decay of purity (1 - linear entropy) for initial cat states in finite Hamiltonian systems 
where an uncoupled part of the Hamiltonian generates regular (integrable) dynamics. 
Such systems are common among theoretical models (e.g. Jaynes-Cummings, ion trap 
quantum computer etc.) and can within certain approximation be even realized in 
the experiments. Our results are relevant also for the decoherence of macroscopic 
superpositions in general systems if the decoherence is faster than any dynamical 
time-scale involved. Note however that a strict decoherence in mathematical sense, 

1. e. irreversible loss of coherence, is possible only in the thermodynamic limit. Still, 
"for all practical purposes" one can have decoherence also in a sufficiently large finite 
system where the actual act of reversal is close to impossible due to high sensitivity to 
perturbations. Almost all theoretical studies of decoherence start from a master equation 
for the reduced density matrix, see e.g. [ü] for the case of harmònic oscillator. Derivation 
of such a master equation is possible only for very simple systems, for instance for a 
harmònic oscillator coupled to an infinite heat bath consisting of harmònic oscillators [ü]. 
Our approach here is different. We do not use master equation but rather start from the 
first principies, i.e. from a Hamiltonian describing system coupled to a finite "bath". 
In the case of regular uncoupled dynamics we are able to derive the decay of purity 
for initial cat states. We stress that the result applies to any integrable dynamics, not 
just e.g. harmònic oscillators, and that the coupling can be quite arbitrary. Thus we 
will describe phenomena that go beyond simple "reversible decoherence" discussed for 
instance for two coupled cavities jjj. Furthermore, our results might shed new light on 
the occurrence of decoherence in infinite-dimensional systems, e.g. harmònic oscillator 
bath P|, by taking the appropriate limits. 

2. Process of decoherence 

Let us write the Hamiltonian of the entire system as 

H = H + 5 - V, H = H c ® l e + l c ® H c , (1) 

where Hq is an uncoupled part of the Hamiltonian and V is the coupling, with 5 being 
its dimensionless strength. We will use subscripts "c" and "e" to denote "central" 
subsystem and "environment" ( "environment" is used just as a label for the part of the 
system we will trace over, without any connotation on its properties). In the present 
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paper we study purity decay for the initial product state of the form 

\m) = ^ (l^cl(O)) + |^(0)>) ® |^ c (0)), (2) 

with all three states |"0 c i,2(O)) and |^ e (0)) being localized Gaussian wave packets. We 
assume that the initial state of the central system is composed of widely separated 
packets (i.e. is a cat state), so that two composing states are nearly orthogonal, 
(^1(0)1^2(0)) ~ 0. However, manipulations with macroscopic superpositions are quite 
nontrivial. In experimental situations the coherence between widely separated packets 
(e.g. | — a) + \a) for large boson parameter a) might be washed out due to presence of 
losses. [Hj- Purity is defined as 

/(*) = tx c p 2 c (t), p c {t) = tr e |V(í))(V(í)|, (3) 

where \i/j(t)) = exp (— iHt/fi)\ip{Q)). For our cat state © the initial reduced density 
matrix reads 

PM = |(|Vci(0)>(^ c i(0)| + |^c 2 (0)>(^c2(0)| + |^i(0))(^ c2 (0)| + |^ c2 (0))^ cl (0)|). (4) 

It has two diagonal terms and two off-diagonal terms also called coherences. These off- 
diagonal matrix elements are characteristic for coherent quantum superpositions. The 
process of decoherence then causes the decay of off-diagonal matrix elements of p c (t), so 
that after some characteristic decoherence time Td ec we end up with the reduced density 
matrix p™ lx (í) which is a statistical mixture of two diagonal terms only, 

Pr X ( r dcc) = -(IV'cl (Tdec)} (V'cl (Tdec) I + | ^c2 faec) ) (^c2 (Tdec) | ) ■ (5) 

The purity of this reduced density matrix is /(rjcc) = 1/2. Whereas the initial density 
matrix p c (0) (/(0) = 1) has no classical interpretation, the decohered density matrix 
p™ lx (Ydec) has. The aim of this paper is to explicitly show decoherence, i.e. to derive 
the transition 

Pc( 0)^p^( rdcc ). (6) 

We are going to do this for systems with a regular uncoupled part H of the 
Hamiltonian. The progress of decoherence will be "monitored" by calculating purity (JHJ) 
which will decay from initial value 1 to 1/2 after t^ cc . When speaking about the decay 
of off-diagonal matrix elements of p c (t) we should be a little careful though as the notion 
of off-diagonal depends on the basis. In has been recently shown (HI CSj that the purity 
in regular systems and for initial localized wave packets decays on a very long time scale 
r p ~ 1/5, which is independent of h. Localized wave packets are therefore very long- 
lived states - the pointer states (they get entangled very slowly), and this is a preferred 
basis which we have in mind when speaking about the off-diagonal matrix elements. 
We will explain the meaning of pointer states more in detail later, after we derive 
analytic expression for purity decay. Provided r p is larger than Td cc we can assume that 
during the process of decoherence propagation of initial constituent states still results in 
approximately product states, i.e. denoting 1^1,2(0) = ex P ( — iiít/^-) |"0 c i,2 (0)) <8> ^(O)), 



Decoherence in regular systems 



4 



we have for t < r p , |^i,2(í)) ~ |"0ci,2(O) ® iV'ei^C*)) - Note that we do not assume 
that \ij) ce (t)) can be obtained by propagation with H c e alone, in fact they can not be, 
because, for instance, the state of the environment l^ei^íO) depends on the state of the 
central system. The above product form of individual states also justifies the form of the 
decohered density matrix p™ lx (jHJ) and is an important ingredient for the self-consistent 
explanation of the process of decoherence (jüj). The condition r p > r dcc means that the 
decoherence time-scale Td ec we are interested in is shorter than the relaxation time-scale 
r p on which individual pointer states relax to their equilibrium. From the result of this 
paper (|T5|) we will see that this condition is satisfied provided the separation between 
two packets is large enough and/or h is sumciently small. 



3. Purity decay 

Let us now calculate off-diagonal matrix elements of p c (t). They are of the form 

tr e |^(í))(^(í)| « \Ai(t))(Mt)\(MtMi(i))- (7) 
Assuming the states of the central system are still approximately orthogonal, 
{^d(t)\Mt)} « 0, the purity is simply I(t) = |(1 + F e (í)), with 

F c (t) = |(^c 2 (í)|^ci(í))| 2 = tr e p el (í)p e2 (í). (8) 

Two reduced density matrices of the environment are p c i,2Íf) = tr c |V ; i,2 (í) ) (^1,2 (í) | - 
Quantity F e (t) gives the size of off-diagonal matrix elements and its decay is an indicator 
of decoherence. It is an overlap on the environment subspace of two states at time 
t obtained by the same evolution from two different initial product states. It is 
similar to the fidelity [TT], which is the overlap of two states obtained from the same 
initial condition under two different evolutions. In fact F c (t) can be connected with 
a quantity called reduced fidelity [12], i.e. the fidelity on the subspace, see also [T3] 
for a connection between decoherence and fidelity. Here we will not use this analogy 
as we will calculate F e (t) directly. Calculating F c (t) is easier in the interaction picture, 
given by |^5f 2 (*)) = M(t)\ip h2 (0)), with M(t) = exp (iH Q t/h) exp (-\Ht/%) being the so- 
called echo operator used extensively in the theory of fidelity decay. Since exp(—iH t/h) 
factorizes, using 2 {t) instead of p c i,2(t) will not change F e (t). For systems where there 
exists an averaging time í ave after which a time average of the coupling in the interaction 
picture, V(t) = exp (iH t/h)V exp (— iHot/Ti), converges, 

V= lim % í T dtV(t), (9) 



T^oo T Jo 

we can further simplify the echo operator. Note that í ave is a classical time if the classical 
limit exists and effective % is sufliciently small, i.e. í ave does not depend neither on 5 
nor on h. It is given by the classical correlation time of the coupling V(t). Nontrivial 
V will typically occur in regular (integrable) systems. But note that only H needs to 
be regular whereas the coupling V (and therefore also H) can be arbitrary. For 5^1 
one can show that the leading order expression (in 5) for the echo operator is 

M(t) = e - i6t? / h . (10) 
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For details see e.g. [Tüj. We proceed with a semiclassical evaluation of F e (t). The average 
V is by construction a function of action variables only and therefore the semiclassical 
evaluation of F e (t) is simplified. The calculation goes along the same line as for the 
evaluation of fidelity [ÜJ and of the purity for individual coherent initial states [TU] . 
A sum over quantum numbers is replaced with an integral over the action space and 
quantum operator V is replaced with its classical limit v(j), where j = (j c ,j e ) is a 
vector of actions having d = d c + d c components if central system and environment 
have d c and d e degrees of freedom, respectively. Let us consider Gaussian packets 
Üa\ipa) oc exp{-(j a -j*)·A o (j o -j*)/2a + i0* • j a /h}, centered at (j*, 0*) where A a is a 
positive squeezing matrix determining its shape. The semiclassical expression of the 
initial density p a (j a ) = (j a |p a (0)|ja) reads 




P«Ü«) = V / detA a exp{-(j Q -j:)·A a (j a -j:)/^}. (11) 

Subscript "a" takes vàlues "cl", "c2" and "e", for three initial packets constituting the 
initial state 0. Writing shortly pi,2(j) = í>ci,2(jc)í>e(je), we have the expression for 
F e (t), 



F e (t) = h- 2d Jdjd]exp 




$ = ü(jc Je) - V(j c , j e ) + u(j c Je) ~ v(j c , j e ). (12) 

Note that the above integral is essentially a classical average (see ^3] for some results) 
over two densities corresponding to two states 1^1,2(0)) and is therefore a sort of cross- 
correlation function. We expand the phase around the position j* of the environmental 
state, $ w [^(j c )-^(j c )](j c -j*)-[^(j c )-^(j c )](jc-j*), where v' e (j c ) = dv(j c , j*)/9j e is 
a vector of partial derivatives with respect to the environment, evaluated at the position 
of the environmental packet j*. Integration over the action variables now gives 

F e (t) = exp (-t 2 /rl c ) , r dec = y/n/C e /ó 

Ce = \ mii) ~ K 1 (v'Mx) - v'M 2 )) ■ (13) 

We can see, that the decay time Td ec for cat state is indeed smaller than the decay time 
for individual states r p ~ 1/6, provided % C e . In the linear response regime C e is just 
the semiclassical expression for 

c c = ^(<[v-<v> c ] 2 > 1+ aF-<F>^^ 

where the subscripts denote with respect to which initial state the average is performed. 
Plugging F e (t) (|T3|) into the expression for purity J(í) we get purity decay for times 
t < r p . We can actually calculate the purity for longer times as well. Namely, after Td ec 
the reduced density matrix is a statistical mixture of states |^ c i(í)) and |^ c2 (í)). The 
purity can then be written as the sum of purities for individual states as there are no 
quantum coherences present anymore. Therefore, a complete formula for purity vàlid 
also for longer times is 

I(t) = \(h(t) + I 2 (t) + 2F c (t)) J (15) 
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where h^if) = tr c p^ 2 (í) are purities of individual states. h^if) decay on a long time 
scale r p and have been calculated in [TQ| . Here we will just list the result 

h,2Ü)= 1 u 1>2 = A^XcK'v'L (16) 

y/det (1 + (ót) 2 u h2 ) 

with a d c x d e dimensional matrix of second derivatives evaluated at the position jj, 
or }*, of the initial state, = 5 2 w(j^ 2 )/9(j c ) fc 9(j e ) l· 

The formula for J(í) (|15|) is our main result. Let us discuss it more in detail. 
First, there are two time scales. A short one on which F e (t) decays and a long one on 
which decay. On a short time scale r dcc purity drops to J(r dcc ) = 1/2, signaling 

decoherence. This decay of I(t) has a Gaussian form l|13j) and is generally faster the 
further apart the centers of two initial states, |^> c i(0)) an d 1^2(0)), are. Expanding 
^(j* 2 ) around the position j* cl we indeed have C c « ±[ü£ • (j c * 2 - j c * 1 )]A~ 1 [< c ' Üc2 ~ &)]■ 
If the coupling v is between all pairs of degrees of freedom we have C e cx d c <i e (j* 2 — j^) 2 . 
The decay time Td cc f)13|) is therefore inversely proportional to the distance between the 
packets |j* 2 — j* x | and to \fd~ c . Similar scaling of decoherence time with the number of 
degrees of freedom has been recently experimentally measured in a NMR context |15j . 
Note that the decoherence time decreases, i.e. we have % <C C e , also when we approach 
the semiclassical limit 7ï — > 0. As the effective h determines the dimensionality of 
the relevant Hilbert space, Td ec will be smaller the larger Hilbert space we have. It 
also decreases with the number of degrees of freedom d e of the environment, meaning 
decoherence gets faster for the environment with more degrees of freedom. Taking the 
thermodynamic limit is not straightforward though. The fact that the decoherence time 
r dcc goes to zero if the number of degrees of freedom d c of the environment is increased 
at the constant coupling 5 is not a surprise. Indeed, if we have a coupling of the same 
strength with infmitely many (and infmitely fast) degrees of freedom, decoherence will 
occur instantly. To have a physical thermodynamical limit one has to take d e — > oo and 
at the same time decrease the coupling strength to higher modes or make the so-called 
"ultraviolet cutoff" as in e.g. Caldeira-Leggett model [Sj. In such a case one would get a 
finite Td cc . For self-consistent description using J(í) (fT5|) one also needs Td cc t avc . Note 
also that decoherence has a Gaussian shape in contrast to master equation approach 
where the decay is exponential. Gaussian form of decoherence has been obtained for 
macroscopic superpositions in [Tïï] . 

Of course, all this holds provided the average perturbation v is nonzero. Interesting 
suppression of decoherence might arise for v = 0, for which extreme stability of quantum 
fidelity has recently been found [TZj. After decoherence time rd cc we are left with a 
statistical mixture of states and further decay of purity I(t), due to the decay of h^it), 
is a consequence of relaxation to final state. This relaxation happens on a long time 
scale r p given by the slower decaying /^(í), see [TÜ] for details. In order to illustrate 
the decay of I(t) we performed numerical simulations. 
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Figure 1. Purity dccay for system 1)17(1 for initial cat state and different í/h — 25, 100 
and 500, full curves from right to left. Dotted curves converging to a singlc curve for 
t > 100 are the theoretical prediction for I(t) 1)15(1 using dccay times 1(19120(1 . Vertical 
arrows indicate times at which we show Wigner functions in Fig. In the inset we 
show the decay for macro-superposition, see text for dctails. 



4. Numerical example 

We consider two coupled anharmonic oscillators, i.e. d c — d e = 1, with the Hamiltonian 

H = j c (ha+a c - A) 2 + ^ e (ha+a e - A) 2 + 6h 2 (a+ + a c ) 2 (a+ + a c f . (17) 

where a+, a a denote boson raising/lowering operators. All initial wave packets are boson 
coherent states, ^^(O)) = e aa+_Q * a |0), where |0) is the ground state. The parameter a 
is chosen as a = with the initial action j* = 0.1 for the environment and j* x = 0.2 

and j* 2 = 0.01 for the two states of the central system. The "squeezing parameter" A 
is for coherent states equal to A = l/(2j*). Other parameters of the Hamiltonian 
are 7 C = 1, j e = 0.6456 and A = 1.2. Coupling strength is set to 5 = 0.01, but note 
that our theory is often not limited to small S ^ü]. Time averaged coupling is easily 
calculated from the classical limit of the Hamiltonian and is v = 4:j c j c . Using this we 
easily evaluate C c (fT3|) 

C e = 16j:(j: i -j: 2 ) 2 ^ 0.058, (18) 

and the decay time Td ec of F e (t) fl~3|) . 

rdec = y/m/m* cl - f c2 \) « 416v^. (19) 
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t=100 



t=300 



Figure 2. (Color onlinc)Wigncr functions for the reduced density matrix p^ 1 (t) and 
initial cat statc Planck constant is % = 1/100 all other data are the same as for 
Fig. ^ On x and y axis we plot \/2h$la and \/2h'^sa, respectively. Disappearance 
of the oscillations is an indicator of decoherence whereas the fading of packets is a 
signature of relaxation. See text for details. 



The matrix u needed in Ii,i{t) (JHJ) is just a number equal to u\p = 64^'^ 2 . Theoretical 
decay of h^if) (H3) is therefore 

7 lj2 (í) = 1 Tpl>2 = r (20) 



In Fig. Qwe show the results of numerical simulation for three different vàlues of h. For 
the smaller vàlues of h we can clearly see the two regimes discussed. Initially J(í) decays 
due to decoherence as described by F e (t). After the decoherence time Td oc , i-e. after 
I(t) falls to the value 1/2, relaxation begins. This regime is described by two decaying 
purities h^if) of individual states and does not depend on h. 

For the largest value of h = 1/25 finite size effects can be observed. The origin 
of finite size effects is twofold: First, the condition of well separated decoherence time 
scale r dcc (fr-dependent) and relaxation time scale r p (/i-independent) results in h <C C e . 
Second, there is a finite saturation value of purity for nonvanishing h, I(t — > oo) ~ h dc . 
For the largest h = 1/25 shown in Fig. the saturation value is / ~ 0.3 which is 
almost as large as purity at the end of decoherence, I(rd cc ) = 1/2. Therefore this 
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finite saturation value will start to influence purity decay already at t^ cc pa 80 and no 
relaxation can be observed. 

In Fig. El we show a Wigner function of the reduced density matrix p^(t) = 
tr c M(í)p(0)M t (t) at times t = 0,50,100 and t = 300. The integral of the square 
of the Wigner function equals purity. At t = we see strong oscillations in the Wigner 
function at the midpoint between two packets. This is a characteristic feature of coherent 
superpositions. At t = 50 we are in the middle of decoherence where the oscillations 
and negative vàlues of the Wigner function have already been reduced. At t = 100 
the process of decoherence has ended, visible as the absence of negative regions in the 
Wigner function. At still longer time t = 300 the relaxation is in full swing as the second 
packet has almost entirely disappeared. In Fig.|2]one can observe that the packets rotate 
around the origin and that the relaxation rates of the two packets are not equal. This 
can be explained by the form of the time averaged Hamiltonian in the interaction (i.e. 
echo) picture, v = 4:j c j e . This Hamiltonian causes rotation of the angle of the central 
system, 9 C = 4j c . Using the value of j* = 0.1 we get that the packets should rotate 
for 90° around the origin in time t ~ 400, which agrees with the data shown. Unequal 
relaxation can be simply understood from the expression for t p i^ (120)1 . The packet with 
smaller j* will relax on a longer time scale. This slower decay is a second order effect 
due to smaller size of the packet in the action direction (i.e. larger A) which causes a 
slower dephasing of angles 6 e and in turn slower decay of purity I(t). 

We should mention that at very large times (i.e. much later after the relaxation) 
we get revivals of purity. These are a simple consequence of having a finite number 
of eigenmodes of the uncoupled system for not too large h, which causes a complex 
beating-like phenomena at large times. Note that these revivals are not simple Rabi 
oscillations between the two subsystems as discussed for instance for cavities in [7] . For 
our data (not shown in the figure) and times t < 10000 we have a revival / ~ 0.95 
at t pa 4000 for k = 1/25, revival of I « 0.5 at t pa 8000 for h = 1/100 and there 
is no revival for H = 1/500 (and t < 5000). The revivals are therefore less prominent 
and happen at larger times for small H. The same would happen if the central system 
would have more than one degree of freedom. For many degrees of freedom systems, 
high sensitivity to perturbations will also effectively prevent the possibility of restoring 
coherence upon time-reversal. 

Finally let us comment on the relevance of our results for the appearance of the 
macro world. It has been pointed out ^B] that the decoherence time for sufficiently 
macroscopic superpositions gets very short. Therefore, the system's dynamics can be 
considered to be regular on this very short time scale, but on the other hand we need 
Tdec 3> ^ave for our theory to apply (í ave does not depend on 5 or Ti). Resolution of 
this problem is immediate if one looks at the quantum expression for C c ()14|) . For 
short times one has to replace V with an "instantaneous" operator V(0) and evaluate 
quantum expectation vàlues using classical averaging. Doing this for our model (fT7j) 
for instance, we get Td ec (í <C ^avc) = Tdcc/4, with Td cc given in Ijlflj) . This theoretical 
prediction agrees with the numerical curve for 5 = 0.06, % = 1/100 and j* 2 = 1, shown 
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in the inset of Fig. ^ where Td oc is much shorter than for the three curves in the main 
plot. Again, the theoretical curve (dotted) agrees with the numèrics (full line). Our 
theory therefore covers regime Td cc 3> í av c as well as rd cc <C í a vc- 

5. Conclusions 

We have derived an analytic semi-classical expression for the purity decay of initial cat 
states in finite composite systems which are integrable in the absence of coupling and 
have a nonzero time averaged coupling. Note that the coupling can in general break 
integrability. Such systems are found both among theoretical models and in experiments. 
Pointer states are identified as localized action packets in the action components for 
which the time averaged coupling in nontrivial. Purity for superpositions of such pointer 
states first decays on a short time scale as a Gaussian, indicating decoherence. On longer 
time scale, after the decoherence, it decays in an algebraic way, signifying relaxation 
to equilibrium. Theoretical results are confirmed by the numerical simulation of two 
coupled nonlinear oscillators. 
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